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XIANPENG HU AND DEHUA WANG 

Abstract. The formation of singularity and breakdown of classical solutions to the 
three-dimensional compressible viscoelasticity and inviscid elasticity are considered. For 
the compressible inviscid elastic fluids, the finite-time formation of singularity in classical 
solutions is proved for certain initial data. For the compressible viscoelastic fluids, a 
criterion in term of the temporal integral of the velocity gradient is obtained for the 
breakdown of smooth solutions. 



1. Introduction 

We are concerned with the formation of singularities in smooth solutions to the multi- 
dimensional partial differential equations of viscoelasticity, especially, of the viscoelastic 
fluids. Viscoelastic fluids exhibit a combination of both fluid and solid characteristics, 
and keep memory of their past deformations due to their "elastic" nature. Viscoelastic 
fluids have a wide range of applications and hence have received a great deal of inter- 
ests. Examples and applications of viscoelastic fluids include from oil, liquid polymers, 
to bioactive fluids, and viscoelastic blood flow past valves; see [14] for more applications. 
For the viscoelastic materials, the competition between the kinetic energy and the internal 
elastic energy through the special transport properties of their respective internal elastic 
variables makes the materials more intractable in understanding their behavior, since any 
distortion of microstructures, patterns or configurations in the dynamical flow will involve 
the deformation tensor. For classical simple fluids, the internal energy can be determined 
solely by the determinant of the deformation tensor; however, the internal energy of com- 
plex fluids carries all the information of the deformation tensor. The interaction between 
the microscopic elastic properties and the macroscopic fluid motions leads to the rich and 
complicated rheological phenomena in viscoelastic fluids, and also causes formidable an- 
alytic and numerical challenges in mathematical analysis. For the viscoelastic materials 
with significant viscosities, the equations of the compressible viscoelastic fluids of Oldroyd 
type ([261 EZ]) in three spatial dimensions take the following form ([5j 151 [15]): 



( 



p t + div(pu) = 0, 

(pu) t + div(/m u) + VP = fiAu + (// + A)Vdivu + div(pFF T ), (1.1) 
F t + u-VF = VuF, 
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where p stands for the density, u £ IR 3 the velocity, and F E M 3x3 (the set of 3 x 3 
matrices) the deformation gradient, and 

P = Apt, A > 0, 7 > 1, 

is the pressure. The viscosity coefficients p, A are two constants satisfying 

p>0, 3X + 2p>0, (1.2) 

which ensure that the operator —pAu — (A + //)Vdivu is a strongly elliptic operator. The 
notation u • VF is understood to be (u • V)F and F T means the transpose matrix of F. 
As usual, we call the first equation in (jl.l|) the continuity equation. For system (|1.1|) . 
the corresponding elastic energy is chosen to be the special form of the Hookean linear 
elasticity: 

W(F) = ±\F\ 2 , 

which, however, does not reduce the essential difficulties for analysis. The methods and 
results of this paper can be applied to more general cases. In the physical regime of 
negligible viscosities, (|1.1|) becomes the system of inviscid compressible elasticity: 

p t + div(pu) = 0, 

(pu) t + div(pu ®u)+VP = div(/)FF T ), (1.3) 
F i + u -VF = VuF. 

We refer the readers to El fT5j [23j [29] for more discussions and physical background on 
viscoelasticity. This paper is devoted to the study of formation of singularities for both 
the inviscid elastic flow (|1.3p and the viscoelastic flow (jl.ip . 

Without the deformation gradient F, the system (jl.ip becomes the compressible Navier- 
Stokes equations. There is a huge literature about solutions to the compressible Navier- 
Stokes equations; see Danchin [6], Feireisl [7], Lions |21| . and references therein. In par- 
ticular, the global weak solutions were constructed in [7J [21] for 7 > |, while in [6] the 
global existence of strong solution was proved in the critical spaces when the solution is a 
small perturbation of the equilibrium. 

When the deformation gradient F does appear, the system (jl.ip is much more compli- 
cated than the Navier-Stokes equations, although the equation satisfied by the deformation 
gradient F is a transport equation and is similar to the continuity equation. Fortunately, 
(jl.ip inhibits a list of local conservation laws which make the analysis for the global 
existence of strong solutions available. The local strong solutions to the compressible 
viscoelastic flow (jl.ip with large data were obtained in Hu-Wang (9[|TT], and the global 
strong solutions to the compressible viscoelastic flow (jl.ip with small data in the Besov 
spaces were established in Hu-Wang |10j 112] and Qiang-Zhang [25] . We remark that for 
these local and global existence of strong solutions, the local conservation laws are cru- 
cial for the dissipation of the deformation gradient, and the property that the curl of the 
deformation gradient is of higher order (Lemma 2.1 in [9]) is also very helpful. For large 
initial data, the global existence of strong or weak solutions for (jl.ip is still an outstanding 
open problem. 

Regarding the global existence of weak solutions to (jl.lj) with large data, among all 
of difficulties, the rapid oscillation of the density and the non-compatibility between the 



FORMATION OF SINGULARITY FOR COMPRESSIBLE VISCOELASTICITY 



3 



quadratic form and the weak convergence are of the main issues; namely, the local conser- 
vation laws are not enough for the convergence of the nonlinear term, especially for the 
term related to the deformation gradient although those terms sit well in the framework of 
div-curl structure, and the higher integrability of the density or the deformation gradient 
is not available up to now. 

For a strong solution to (jl.ip with large data, although the local existence was proved 
in [9], we expect that it will break down in a finite time as for the strong solution to the 
compressible Navier-Stokes equations ([1] [28l 113]). The breakdown of smooth solutions 
is due to the lack of control of the hydrodynamic variables, for instance, the L°° norm 
of the gradient of the velocity or the L°° norm of the density. As it is well-known, the 
L°° norm of the gradient of the velocity controls the L°° norm of the density and the 
deformation gradient in the compressible viscoelastic fluids (jl.ip . In this paper, we will 
provide a criterion and justify the blowup phenomena. Comparing with the compressible 
Navier-Stokes equations, the viscoelastic flow (jl.ip is more complicated and thus more 
delicate and new estimates are needed for the analysis of strong solutions. More precisely, 
the main difficulty lies in the estimates of the gradients of the density and the deformation 
gradient, and the estimates on the L^H^. bounds of Vp and VF are crucial. We find a 
criterion for breakdown of strong solutions of (jl.ip in term of the temporal integral of the 
L°° norm of the velocity gradient. 

For the inviscid flow (jl.3p . similar to the compressible Euler equations (Sideris [31J), we 
expect that the smooth solution to the system (jl.3p will develop singularities in a finite 
time. We will first reformulate the system (jl.3p into a symmetric hyperbolic system so 
that a local smooth solution can be obtained from |16j I24j. Then we will prove that the 
smooth solution cannot exist globally in time under some restrictions on the initial data; 
that is, the finite-time formation of singularities is essentially inevitable provided that the 
initial velocity in some region near the origin is supersonic relative to the sound speed at 
infinity. The proof will follow the idea of Sideris |31] with more subtle estimates on the 
deformation gradients. 

For the incompressible viscoelastic flows and related models, there are many papers 
in literature on classical solutions (cf. [21 El [TTl [TH1 ED] and the references therein). On 
the other hand, the global existence of weak solutions to the incompressible viscoelastic 
flows with large initial data is also an outstanding open question, although there are some 
progress in that direction ( [TH [22], [23] ) . For the inviscid elastodynamics, see [32] and their 
references on the global existence of classical solutions. 

The rest of this paper is organized as follows. In Section 2, we will explain the mechanism 
which will ensure a local existence of smooth solution to the inviscid compressible elastic 
fluid and also provide a proof of the finite-time formation of singularities. In Section 3, 
we will consider the compressible viscoelastic fluids and prove a blowup criterion in term 
of the the L°° norm of the gradient of the velocity. 

2. The Inviscid Case 

In this section, we consider the formation of singularities in smooth solutions of the 
inviscid flow (11.3P in M 3 with sufficiently smooth initial data: 

(p,u,F)\ t=0 = (po{x),u (x),F (x)), i£R 3 . (2.1) 
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We assume that 

po(x) > for all x £ M 3 , 
there exist positive constants po and R > such that 

(po(a:),uo(a;),Fo(a;)) = (p ,0, /) for all |x| > R, 
where I is the 3x3 identity matrix, and 

dw(p Fl) = 0, 



(2.2) 
(2.3) 
(2.4) 



One useful property of the deformation gradient F is the following (see Lemma 6.1 in 
101): 



Lemma 2.1. If (p, u, F) is a smooth solution of (|1.3p . and p,F initially satisfy ()2.4|) . then 
the following identity holds for any time: 

div(/)F T ) = 0. (2.5) 

Under the assumption (|2.4p and using Lemma |2. 11 we can rewrite the system fjl .3j) for 
smooth solutions (p, u, F) with p > as: 

1 dP . , 
- f - + d 1 vu = 0, 



dFi 
I dt 



+ VP = F ;7 ,V, ,F, 7 , := ^2(F k • VF A 
= VuF, 



(2.6) 



k=l 



where 



d d „ 



c 2 =J 4 7(7 _ 1)p 7-2 ) p 

and F,; is the i-th column of the matrix F. Set 
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1,2,3, 



where {e\, e 2 , e%} is the standard basis of R 3 , I12 is the 12 x 12 identity matrix, ej is the 
transpose of e^, and / is again the 3x3 identity matrix. Then, in view of (12. 6p . the system 
(jl.3p can be written as a symmetric hyperbolic system of the form 



(2.7) 



4 = 1 
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with Aq(V) > and Ai(V) is a 13 x 13 symmetric matrix for each i = 1,2, 3. According to 
the well-known result in \16\ 124] . the hyperbolic system of conservation laws (|2.7p admits 
a local C 1 solution on some time interval [0, T), provided the initial data are sufficiently 
regular; moreover, p > on R 3 x [0, T). 

Since po(x) = po, uq(x) = 0, and Fo(x) = / for all \x\ > R, then a, the sound speed at 
infinity, is given by 

The following proposition is an immediate consequence of local energy estimates (cf . {30] ) . 
It simply states that the maximum speed of propagation of the front of a smooth distur- 
bances is governed by a. 

Proposition 2.1. I/(p,u,F) G C^M 3 x [0,T)) is a solution of flL3]) and (j27l|) . then 
(p, u, F) = (p , 0, 1) for all \x\ > at + R and < t < T. 

As for the compressible Euler equations (Sideris [31]), we expect that the smooth so- 
lution to the system (jl.3p will develop singularities in a finite time. The result in the 
following theorem shows that the C 1 solution to (jl.3p and (|2.ip does not exist globally in 
time under some restrictions on the initial data; more precisely, it states that the finite- 
time formation of singularities in the three-dimensional inviscid compressible elastic fluid 
is essentially inevitable provided that the initial flow velocity, in some region near the 
origin, is supersonic relative to the sound speed at infinity. 

In order to state our result, we define 

m (t) = / (pOM) - Po)dx, 
Jm.3 



•F{t) = / p(x, t)x • u(x, t)dx, 
Jrz \ 2 2 7 - 1 / 

with 

P = Pifio) = Apl 

and 

D{t) = {x G M 3 : |x| < at + R}. 

Remark 2.1. Proposition 12.11 implies that the integrands of m(t), J~(t), and £(t) are 
identically equal to zero outside D(t). 

For two 3x3 matrices A and B, the following notations will be used: 

3 3 

A : B = AijBij, \A\ 2 = A\y 

i,j=l i,j=l 

Lemma 2.2. Let (p, u, F) be a C 1 solution of (jl.3p with initial data ()2.ip - (j2.4p . Then 
£{t) is conserved; that is, 

S'(t) = 0, £(t)=£(0), (2.8) 

for all t > 0. 
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Proof. Multiplying the second equation in (jl.3p by u and using the first equation in (jl.3 
yield 

d r (\ , |2 p-Pq 



, , , -plul 2 + )dx = - pFF ' : Vudx. (2.9) 

at Jr3 \2 7 - 1 J J RS 

On the other hand, from the third and then the first equations of (jl.3p . one deduces that 

= ir\ F \ 2 + 2 P F : ( VuF - u • VF ) 



dV 

F| 2 + 2pF : (VuF) - pu • V|F| 



at 

^|F| 2 + 2/9F : (VuF) + div(/m)|F| 2 - div(pu|F| 2 ) 
2pF : (VuF) - div(pu|F| 2 ). 



Since 



plF-Il^p^. + p^-l) 2 

i^j 1=1 

_ 3 
= ^F^-2^F !1 +3/ ) 

i,j i=l 

= p\F\ 2 -2ptr(F-I) -3p, 
where tr(F — I) denotes the trace of the matrix F — I, then 

HF| 2 = p\F - I\ 2 + 2ptr(F - /) + 3(p - p ) + 3p , 

thus 

^ (p\F - I\ 2 + 2ptr(F - /) + 3(p - po)) = 2pF : (VuF) - div( / ou|F| 2 ). 

Integrating the above equality, we arrive at 
1 d f* 

-- / (p\F - I\ 2 + 2ptr(F - /) + 3(p - po)) dx 
z at j R 3 

= pF : (VuF)dx = / pFF T : Vudx. 
Adding (|2.9|) and (|2.10p together, one has 

Jt L + 5* - /|2 + " tr(F - " + 1 ( " - *> + ^t) dl = °' (2J1) 

We note from the first equation of (j 1 . 3 [) that 

m '(t) = A / (p-po)dx = - [ div(pu)dx = 0, (2.12) 

thus, 

m(t) = m(0). (2.13) 



(2.10) 
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Due to Lemma |2.H we have, using integration by parts, 

pF T : Vudx = 0. 

From the first and third equations in (|1.3p . it is easy to deduce that 
<9 t (ptr(F - I)) + div(ptr(F - I) u) = pF T : Vu. 
Integrating the above equality yields 

/ ptr(F - I) dx = 0, (2.14) 



d_ 

dt 



and then 



/ pti{¥ - I)dx = i po tr(F - /) dx. (2.15) 
Substituting (|2~12]) and (pEUD into (l2~TT]l . we have 



/ /I , ,2 , 1 lr r|2 , P-Pq 



-p|u| 2 + -p|F - I\ 2 + ^ dx = 0. (2.16) 

dt J-^3 \z 2 7-1 J 

Therefore, (I2.8P follows from (|2.16p . The proof is complete. □ 
Now we state and prove the finite-time formation of singularity for (jl.3j) . 



Theorem 2.1. lei (p, u, F) G C^M 3 x [0, T)) 6e a solution of (fL3l) wii/i zmiia/ data 

m(0) > 0, (2.17) 

J-(O) > ^ai? 4 ||p |k-, (2.18) 

and 

/ po tr(I - F )dx > 2£(0), (2.19) 

then T is necessary finite. 

Proof. Let (p,u,F) be a C 1 solution of and J3U. Then, J 7 G C^O.T) and 

J r '(t)= / x • (ptu + put)dx. 

It follows from the Proposition 12. 1( the first two equations in (II. 3p . and integration by 
parts, that 

F'{t)= p\u\ 2 dx + 3 (P-P )dx + x ■ div (pFF T ) dx. (2.20) 

JD(t) JD(t) JD(t) V 7 
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Using the Holder inequality, f|2. 13|) and (|2,17p . we have, 



Pdx = A / p~ f dx 

D(t) JD(t) 



> A(vQLD(t)) 1-7 ( [ pdx] 

\JD(t) J 

= ^(voLD^)) 1 " 7 (m(0) + voLD(t)A) 

> vo\D(i)Apl = / P dx, 

JD(t) 



thus 



/ (P - P )dx > 0. (2.21) 

JD(t) 

Note that Lemma |2. II implies 

div(pFF T ) = div(/)(F - J)(F - /) T ) + div(p(F - /)). 
Then the divergence theorem yields, using Lemma 12.21 (|2.15|) . (|2.19|) and (|2.2ip . 



(2.22) 



/ x ■ div(pFF T )dx 
JD{t) 

= / x ■ div(/>(F - J)(F - I) T )dx + / x ■ div(p(F - I))dx 

J D{t) JD(t) 

= - [ I : p(F — I)(F - I) T dx - [ I : p(F - I)dx 

J D{t) JD{t) 

= - p\F-I\ 2 dx- / ptr(F-I)dx 

JD(t) JD(t) 

> -2£{t) - / ptr(F- I)dx 

JD(t) 

= -25(0) + Po tr(I - F )dx > 0. 
Therefore, I^M)-^22^ yield 

J"'(f) > / p\u\ 2 dx. (2.23) 

JD(t) 
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On the other hand, (|2,23p and the Cauchy-Schwarz inequality lead to 



where 



> F{t) 2 (at + R)- 2 if pdx) 

\JD(t) J 

= F{t) 2 (at + R)' 2 I [ p (x)dx) 
\JD(t) J 

> ^ (f7t + jR )5|| po || ioo ^ jr (t) 2 j 

IIa)IU°° = su p Poi x )- 



Since J"(0) > 0, then (HEMl) implies that F(t) > for < t < T, and that 
^(O)- 1 > .F(O)- 1 - F{T)- 1 

1 1 



(2.24) 



- 167TCr||po||L°° V^ 4 (o-T + i?) 4 / 

which shows that T cannot become arbitrarily large without contradicting the assumption 
(I2,18p . Therefore the proof is complete. □ 

Remark 2.2. Taking po(x) = po, we see that (|2.1T|) - f)2.19[) can be satisfied if uo is super- 
sonic relative to the sound speed a, the diagonal entries [/ — Fq]h G (0, 1), i = 1, 2, 3, and 
po is not very large. 



3. The Viscous Case 

This section is devoted to the study of formation of singularity and breakdown, espe- 
cially of the blowup criteria for the smooth solutions of compressible viscoelastic fluids 
(jl.ip in R 3 with sufficiently smooth initial data: 

(p,u,F)\ t=0 = (p (x),u o (x),F (x)), xeR 3 . (3.1) 

The goal is to obtain a blowup criteria in term of the L°° norm of the velocity. In other 
words, the L°° norm of the velocity controls the blowup for the compressible viscoelastic 
fluids. For this purpose, we need to assume that 

7p > A. (3.2) 

Obviously, this condition will be fulfilled physically if viscosities p, and A satisfy the con- 
dition (jl.2p and A < simultaneously. 
Denote D k and D§ as: 

D k (R 3 ) := {/ G Lj oc (R 3 ) : ||V fc /IU W < oo}; 

and 

D k (R 3 ) := {/ G L 6 (R 3 ) : ||V*/[|^(R3) < oo}. 
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Now, our blowup result for the system (jl.ip takes the following form: 
Theorem 3.1. Assume that the initial data satisfy 

< po G # 3 (K 3 ), u G Dq (M 3 ) n Z) 3 (R 3 ), F G F 3 (M 3 ), div(poFj) = (3.3) 

and 

- ^Auo - (A + /i)Vdivu + AV/$ = p 5 (3.4) 

/or some g G i? 1 (M 3 ) with ^fp^g G L 2 (M 3 ). Lei (p, u, F) be a classical solutions to the 
system (jl.ip satisfying 

r( P ,F)eC([o,T*],H 3 (R 3 )), 

u G C([0, T*], Dj(l 3 ) n D 3 (R 3 )) n L 2 (0, T* ; D 4 (S 3 )), 
u t G L oo (0,T*;^(R 3 ))nL 2 (0,T*;D 2 (M 3 )), 
{y/put GL°°(0,T*;L 2 (R 3 )), 

and T* 6e t/ie maximal existence time. If T* < oo and (13. 2p holds, then 



(3.5) 



lim / || Vu|| L ooC] K 3)di = oo. (3-6) 

T->T* Jq 

The local existence up to a possible finite time T* can be constructed in the functional 
framework (|3.5p with the compatibility condition (|3.4p in spirit of the corresponding results 
for the compressible Navier-Stokes equations (see for example [3]), thus we omit the details 
in this paper and focus on the breakdown of smooth solutions. To prove Theorem 13. 1\ the 
main difficulty lies in the estimates of the gradients of the density and the deformation 
gradient provided that the quantity in (]3.6p is finite. In fact, the key estimates in our 
analysis is L^H^. bounds of Vp and VF. 

3.1. Regularity of solutions. In this subsection, we will derive a number of regularities 
of solutions (p, u, F) to the system (II. lh provided 

f T 

lim / || Vu | |^oo(R3 )dt < oo. (3-7) 

T—tT* Jq 

The standard energy estimate gives 

sup f||Vpu|| 2 2 + ||p|| L ^ + \\^pF\\ 2 L2 ) + [ T \\Vu\\ 2 L2 dt <C, < T < T*. 

0<t<T v 7 Jo 

Moreover, from the two transport equations for the density and the deformation gradient, 
the assumption (|3.T[) yields the L°° bounds for the density and the deformation gradient 
respectively. 

The first step for the regularity of solutions is to improve the integrability of the velocity. 
In fact, we have 



Lemma 3.1. Under the assumption (13. 2p . there exists a small 5 > such that 

sup / p\u\ 3+5 dx <C, < T < T*, 

o<t<T Jm.3 



where C is a positive constant depending only on ||p||l°° an d ||F||l°°- 
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Proof. The argument is similar to that of |13j . For q > 3, we multiply the second equation 
in (jl.ip by g|u| IJ ~ 2 u and integrate over M 3 to obtain, using the conservation of mass, 

p\u\ q dx+ [ fg|u|"- 2 (/x|Vu| 2 + (A + ^)|divu| 2 + /i(g-2)|V|u|| 2 ) 

at J R 3 JR3 \ 

+ q(X + /i)(V|u| 9 ~ 2 ) • udivujdx 
= Aq I div(|u|^ 2 u)p 7 dx -q I pFF T : V(\u\ q - 2 u)dx ,„ on 

Jr 3 Jr 3 

<C f ^p\u\ q - 2 \Vu\dx < e f \u\ q - 2 \Vu\ 2 dx + C [ p\u\ q - 2 dx 
Jr 3 Jr 3 Jr 3 

9-2 

<e [ \u\ q - 2 \Vu\ 2 dx + C ( [ p\u\ q dx] " , 

Jr 3 \Jr 3 J 

for e > 0. As in [13], if 7p > A, one has 

g|u| 9 ~ 2 (V|Vu| 2 + (A + /i)|divu| 2 + p(q - 2)|V|u|| 2 



+ q(\ + p){X7\u\ q - 2 ) -udivu (3-9) 
> C|u| 9 - 2 |Vu| 2 . 

Substituting (|3.9p into (|3.8p . and taking e small enough, we get by GronwaM's inequality, 

sup / p\u\ 3+6 dx <C, < T < T* , 

o<t<T Jr 3 

for some small 5 > 0. □ 

The second step for the regularity considers the integrability of the material derivative 
of the velocity, which is useful for the bounds on Vp. To this end, we have, 

Lemma 3.2. Let 

G = pvn + pu ■ Vu. 

Then, 

[ [ G 2 dxdt+ sup / \Vu\ 2 dx<C [ [ (|Vp| 2 + \VF\ 2 )dxdt + C 

JO Jr 3 0<t<TjR3 J Q J R 3 

for allO<T < T* . 

Proof. From the L°° bound of p, we have 

r-T r r-T r- r-T 



nG 2 dxdt<c[ I pu 2 dxdt + 2 [ [ \pu • Vu\ 2 dxdt. 
_: 3 Jo Jr 3 Jo Jr 3 



(3.10) 



12 XIANPENG HU AND DEHUA WANG 

We proceed to estimate the right-hand side of (|3.10p term by term as follows. For the last 
term of (|3.10p . we have 



T 

JR 3 



p 2 \u\ 2 \Vu\ 2 dxdt 



<C [ T \\^pu\\l 4 \\X7u\\ 2 L4 dt 
Jo 

<C f \\^u\\l 3+ 4^u\fe a \\Vu\\ 2 L4 dt 
Jo 

r T 

< C / l|Vulfc a l|Vul|r<»<ft 



(3.11) 



o 



T 



where 



< C sup ||Vu||^2 Q / \\Vu\\ L °odt 

0<T<T* JO 

< C sup ||Vu||^2 Q , 

0<T<T* 



a 2-a 1 , 3+5 „ 
-, 1 < a = < 2. 



3 + 5 6 2' 3-6 

For the first term, we multiply the second equation of (jl.ip by ut and integrate over 
to obtain 



2 , f _ . d f ///._ l2 A + /X ,o\ . 

pu t dx + / pu • vu • u t dx + — / I— |Vu| H — |divu| lax 

Jr3 at J R 3 V 2 2 / 

= A p^divUfdx — / /9FF T : Vu t dx. 

Jr 3 Jr. 3 

From the first equation and the third equation in (jl.ip . we deduce that 



(3.12) 



p^divutdx 

/ //Mivudx — 7 / ptp 7_1 divudx 
= — / p 7 divudx + 7 / div(pu)/9 7_1 divudx 

= — I p 7 divudx- / p 7 u • Vdivucte + (7 - 1) / p 7 |divu| 2 dx ( 3 ' 13 ) 

/ K 3 J K 3 

d f 

< -j t y^ 3 P 7 divudx + C||^u|| L2 ||V 2 u|| L 2 (IR 3 ) + C||Vu||£ 2(R3 ) 

-JtL pMWudx + c + £ \\ g W 2 l'^ + c W v p\\hm 

+ C\\VF\\ 2 +C\\Vu\\ 2 



d_ 

dt 
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and 



/ pFF T : Vu t dx 

Jm. 3 

= — [ pFF T : Vudx - [ d t {pFF T ) : Vudx 



< — pFF : Vudx — / pFikFj^u ■ VdjU l dx 

dt Jr3 J R 3 

+ C f p|F| 2 |Vu| 2 dx 

< - / pFF T : Vudx + C\\ Vpu|| i2 ||V 2 u|| L2(R 3) + C||Vu|| 2 2 
dt J R 3 

-JtJ 3 pFpT : Vudx + c + e||G|1 + c W v p\\hm 
+ c|I vf IIl 2 (r 3 ) + c|I Vu IIl 2 (r 3 )' 

where we used 

d t (pFF T ) + u • V(pFF T ) = Vu(pFF T ) + pFF T (Vu) T - pFF T divu. 

Note that 



(3.14) 



JO Jm. 3 



pu ■ Vu • Utdxdt 



< 



npu 2 dxdt + I i p\u ■ Vu\ 2 dxdt 
jjs Jo Jrs 



< 




3-a 



pu t dxdt + C sup ||Vu||r2 
lo Jr 3 0<T<T* 



Substituting (13.131) and (13.140 to (|3.12j) . using Gronwall's inequality and the Cauchy- 
Schwarz inequality, we have, 



sup ||Vu|| l2(R 3) < C; 

0<t<T* 



and 



pu t dxdt 

rT 



<C [ [ (\Vp\ 2 + \VF\ 2 )dxdt + s [ [ G 2 dxdt ( 3 - 15 ) 
Jo Jr 3 Jo Jr 3 



+ C sup ||Vu|| 3 .2 Q + C. 

0<T<T* 

Here we used the following estimates: for all < t < T 



and 



pFF : Vudx 
p^divudx 



<C7||VpF|| l2 ||Vu|| l2 <C7 + ^||Vu|| 2 2 



^CVIMlVuH^ <C + £||Vu||£ 2 . 
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Choosing a sufficiently small e, from (|3.10p and (|3.15p . one has 

f [ G 2 dxdt+ sup f \Vu\ 2 dx<C [ [ (|Vp| 2 + \V¥\ 2 )dxdt + C. 

JO JR 3 0<t<TJR3 J J R 3 

The proof is complete. □ 

With the aid of the estimate of the material derivative of the velocity, we can obtain 
the L°°(0, T; L 2 ) estimates of Vp and VF and we actually have 



Lemma 3.3. Under the assumption (|3.2p . the following estimates hold for < T < T* : 

sup / \Vp\ 2 dx < C, (3.16) 

0<T<T* Jm.3 



sup / |VF| 2 dx < C, (3.17) 

0<T<T* 



npu 2 dxdt + sup / |Vu| 2 cte + / ||u||^2( R 3\c£i < C, 
_;3 o<t<t* Jm? Jo v ' 



(3.18) 



Proof. The arguments of (|3.16p . (|3.18p are similar to Proposition 2.4 in [13] provided 
(I3.17P holds, and thus we omit them. For (|3.17p . we can proceed as follows. 

Differentiating the third equation in (jl.ip with respect to X{ (denoting d Xi by d{) and 
multiplying the resulting identity by 2<9jF yield 

<9t|<9;F| 2 + u • V|diF| 2 = 2V5juF : <9;F + 2Vu<9;F : <%F. 

Integrating the above identity over M 3 , one has 

d 



dt 



diF\ 2 dx = / |<9jF| 2 divudx + 2 1 (V^uF : <%F + Vud;F : 0»F) o?x 

< C||Vu|| LO c||aiF||| 2 +C||aiF|| L 2||V 2 u|| L 2 ( 319 ) 

< C||Vu|Uoo||c>,F||| 2 + C||^F|| 2 2 + C||V 2 u|| 2 2 

< ciivuHlooIi^fii 2 , + c\\di¥\\ 2 L2 + c\\d lP \\ 2 L2 + c||G|| 2 2 



since 



-/iAu - (A + /i)Vdivu = —G — AVp" 1 - div(pFF T ). 
From (13.16p . (I3.19P and Gronwall's inequality, we obtain (|3.17p . The proof is complete. □ 

The next step is to obtain the uniform in time estimate for as the following: 
Lemma 3.4. Under the assumption (|3.2p . the following estimates hold for allO <T < T* , 



sup \\^/pu t \\ 2 L2 + [ [ \Vu t \ 2 dxdt < C, (3.20) 

0<t<T Jo Jk. 3 

sup ||u|| H 2 < C. (3.21) 

0<T<T* 
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Proof. The argument of (|3.2ip is similar to Proposition 2.5 in [13], and thus we omit it. 
To prove (|3.20p . we differentiating the second equation in (jl.ip with respect to t to obtain 



putt + pu • Vu t — /iAii( — (A + /i)Vdivu( + AV(p" 1 
= -p t (u t + u • Vu) - pu t • Vu + div(pFF T ) t . 



(3.22) 



Taking the inner product of the above equation with in L 2 (R 3 ) and integrating by 
parts, one obtains, 

I— I pufdx+ I (/x|Vui| 2 + (A + ^)|divu t | 2 )dx 
2 at j R 3 J R 3 

-A I (p^tdivutdx + / (pFF T ) t : Vu ( (3.23) 

Jr 3 Jr 3 

pu ■ V (^ u t + u • Vuj u t + pu t ■ Vu • u t J dx. 
From the first and the third equations in (jl.ip . we have, 



since 



and 



(p 7 )tdivu t dx 

= — I div(p 7 u)divuf(i3; — (7 — 1) / //Mivudivufdx 

JR 3 ' J R 3 

= / (p 7 divu + V/9 7 • u)divu 4 dx — (7 — 1) / //Mivudivu^x 

Jr 3 Jr 3 

< C||Vu t || L 2||Vu|| z2 +C||u|| i -||V/9|| L 2||Vu t || i 2 

< e||Vu t ||| 2 + C||u|| 2 ?2 

(p 7 ) 4 + div(p 7 u) + (7 - l)p 7 divu = 0; 



(pFF T ) f : Vu t dx 

3 

- / u • V(pFF T ) : Vu t dx + / VupFF T : Vu t dx 

Jr 3 Jr 3 

+ / pFF T (Vu) T : Vu t dx + / pFF T divu : Vu 

Jr 3 Jr 3 



t dx 

< C||u|| L o (||Vp|| x2 + ||VF|| L2 )||Vu t || L2 +C||Vu|| L 2||Vu t || L 2 

< e||Vu t || 2 2 +C\\u\\ 2 H2 . 



(3.24) 



(3.25) 
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For the right-hand side of (|3.23p . we have 



pu • V 



-Ut + u- Vu|u( 



+ pu t • Vu • ut ) dx 



< 



/ fp|u||ut||Vut| +p|u||ut||Vu| 2 + p|u| 2 |ut||V 2 u| 
+ p|u| 2 |Vu||Vut| + p|u t | 2 |Vu|W 



< C(J|u|| L 6|| v / put|| L3 ||Vut|| L 2 + ||u|| L 6||ut|| L 6||Vu|| L 3 
+ ||u 2 || L 3||ut|| L 6||V 2 u|| L 2 + ||Vut|| L 2||Vu|| i 6||u 2 || L 3 



+ H^/putll^allVuHio 

< C(j|Vput||| 2 ||Vut||2 2 + ||Vut|| L2 ||Vu|| L 6 + ||ut|| i 6||V 2 u|| L2 
+ ||Vu|| L 6||Vut|| L 2 + || v ^pu t ||| 3 ||Vu||i / o J 

< e||Vut|| 2 2 + C(l + ||Vu|| L oo)||^/put|| 2 2 + C||u||^2. 
Substituting SSTHty . (13351) and (EOHD into ([323]), we get 



1 



-pu t dx + p / |Vut| dx 



d 
dt 

< e f \Vu t \ 2 dx + C((l + ||Vu|| L ^)||^ut|||2 + Hull^). 

Thanks to the compatibility condition, it holds 

^u t (0,x) G L 2 (R 3 ), 
and thus, for sufficiently small e, (|3.27p gives us (|3,20p . The proof is complete. 



(3.26) 



(3.27) 



□ 



Finally, we now can obtain bounds of the first order derivatives of the density and the 
second derivatives of the velocity. 



Lemma 3.5. Under the assumption (|3.2p . the following estimates hold for < T < T* , 

(3.28) 



sup (||pt(i)|| L 6 + IIpII^m) < C, 

0<t<T 

sup (\\F t (t)\\ L e + \\F\\ wh6 ) <C, 
0<t<T 



U 



w 2 - 6 



dt < C. 



(3.29) 



(3.30) 



Proof. The arguments of (|3.28p and (|3.30p are similar to those in [13] provided (|3.29p is 
verified, and hence are omitted here. To establish (|3.29p . we first deduce from the previous 
estimates that, 

ut € L 2 (0,T;L 6 (M 3 )), G G L 2 (0, T; L 6 (M 3 )). 



FORMATION OF SINGULARITY FOR COMPRESSIBLE VISCOELASTICITY 17 

Differentiating the third equation in (jl.ip with respect to Xi, and multiplying the resulting 
identity by 6|<9jF| 4 <9jF, one has, 

= -6 / \diF\ 6 divudx - 6 / \diF\ 4 diFdiVuFdx - 6 / \df\ 4 VudiFdiFdx ( 3 - 31 ) 

JR3 J R 3 J R 3 

< c(||Vu|| L oo||VF||| 6 + ||VF||| 6 (||Vp|| L 6 + ||VF|| L 6 + ||G|| £8 ; 



The above inequality and GronwalFs inequality imply 

sup ||VF|| L 6 < c. 

0<t<T 

This estimate, together with the third equation in give the bound of Ft in (|3,29|) . 

The proof is complete. □ 

3.2. Proof of Theorem l3.ll To begin with, we have the following higher order estimates 
for the density and the deformation gradient: 

Lemma 3.6. For all <T <T* , the following estimates hold: 

\\p\\l°°(h 2 ) + \\Pt\\L°°(m) + IIpmIIl 2 < C, 

IIp 7 IIl-»(^2) + ll(p 7 )illL-»(Hi) + II(p 7 )«IIl2 < c, 

and 

\\F\\l°°(h 2 ) + llFtlli 00 ^ 1 ) + II f «IIl 2 < c. 

Proof. The argument for p and p 7 is similar to those in |13j . and thus we will focus on the 
estimates on F. For this purpose, we first observe that from the second equation in (jl.ip . 
we have 



Nina < C(\\G\\ m + ||V^|| H i + ||div(pFF T )|| H1 ^ 

< c(||G|| ff i + HvVIIl 2 + C + \\VF\\i 4 + [|Vp||i4 + ||V 2 F|| L2 ) 

< c(\\G\\ H i + HvVIIl 2 + C + ||V 2 F|| i2 X 



since 



and 



1 3 



\\Vp\\t*<\\Vp\\l 2 \\Vp\\l 6 <C 



|VF||| 4 < ||VF||£ 2 ||VF||* 6 < C. 



Applying V to the third equation of (jl.ip to yield 

(VijF)t + u • V(V fj F) + V;u • VVjF + Viju ■ VF 
= (V ii Vu)F + VuVyF + VjVuVjF + V^VuViF. 
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Multiplying the above identity by 2V«F and integrating over M 3 , one obtains 



i / |V 2 F| 2 <& 



dt . 

< C / f |V 3 u||F||V 2 F| + |Vu||V 2 F| 2 + |V^u||VF||V^F| )dx 



+ I u-V|V 2 F| 2 dx 



< C f f|V 3 u||F||V 2 F| + |Vu||V 2 F| 2 + |V 2 u||VF||V 2 F|W 

Jr* v 1 

- / divu|V 2 F| 2 cte 
<C I ( |V 3 u||F||V 2 F| + |Vu||V 2 F| 2 + |V 2 u||VF||V 2 F|W 



< C( ||Vu|M|V 2 F|| 2 2 + ||V 2 F|| 2 2 + \\G\\ 2 H1 + ||VV|| 2 2 + C+ ||V 2 u 112 



L6 

2t? 1 1 2 i s~i / 'II /"f 1 1 2 I || V72,, 1 1 2 



< C(||Vu|| L oo + 1)||V 2 F||£ 2 + C(||G||^i + \\V 2 u\\ L * + 1) 



since 



||VF|| L a < ||VF||| 2 ||VF||J 6 <C. 
Because u E L 2 W^' 6 , by Gronwall's inequality, the above inequality yields the bound 



< c. 



The estimates of F^ and F^ follow from the similar energy methods and the third equation 
of (jl.ip . The proof is complete. □ 



When we state the H 3 regularity of the solution (p, u, F), the following estimate is 
useful. 



Lemma 3.7. For < T < T* , we have 



/ / pu 2 t dxdt + sup / \Vu t \ 2 dx < C. 

JO Jm? 0<t<TjR3 
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Proof. Differentiating the second equation in (jl.ip with respect to t, muitipiying the 
suiting equation by Uu, and then integrating over M 3 , one obtains 



pu\ t dx + / pu ■ Vu ( • u tt dx + 



A / (p^tdivuttdx - / (pFF T ) t :V u tt dx - / p t {u t + u • Vu)u tt 

JR3 JR3 J R 3 



divu+| 2 ) dx 



dx 



dx 



dx 



pu t • Vu • uttdx 

= A^- [ (pi)tdivutdx -A I (p") tt divu t dx - ^ / (pFF T ) t : Vu f 

+ / (pFF T ) tt :Vu t dx - [ pu t ■ Vu • u tt dx - — [ -p t \u t 
Jr 3 j R 3 at J R 3 I 

+ 77 / PttWt\ 2 dx - — [ p t (u ■ X7u)u t dx + [ pttVL • Vu • u t dx 

1 Jr3 at Jm.3 J R 3 

+ / p t u t ■ Vu • u t dx + p t u- Vu ( • u t dx. 

Jm 3 Jm 3 

Observe that 



pu • Vu ( • uttdx 



< e\\^u t t\\ 2 L 2 + C||Vu t ||f 2 ; 



pu t • Vu • u tt dx 



< e\\^u t t\\ 2 L 2 +C||Vpu t ||| 3 ||Vu||| 6 

< e\\^/pu t t\\ 2 L 2 + C\\^/put\\ L 2\\u t \\ L 6 
<£\\^putt\\ 2 L 2+C\\Vut\\ L 2- 



ptt\ut\ dx 



/ div(p t u + pu t )\u t \ 2 dx 

Jm, 3 



since 



implies that 



= / (p t u + pu t ) ■ V\u t \ 2 dx 

Jm. 3 

< C\\pt\\ L 6\\u\\ L (i\\u t \\ L (i\\Vut\\ L 2 + / p\u t \ 2 \Vu t \dx 

Jm 3 

< C||Vut||| 2 + C , || v /put|| L 2|| x /pu t || L 6||Vu t || i 3 
<C||Vu t ||| 2 +C||Vu t || L2 ||u t || H2 

< C||Vu t ||2 2 + C\\Vu t \\ L 2(\\^putt\\ L 2 + ||Vu 4 || i2 + 1) 
<C\\Vu t \\ 2 L 2 + e\\^utt\\ 2 L 2+C 

pAu t + (A + p)Vdwu t = G t + AV{pi) t - div(pFF T ) t 

||u t || ff2 < C(\\G t \\ L 2 + \\V(p^)t\\ L 2 + ||div(pFF T ) t || i2 ) 

< C(\\^pu t t\\ L 2 + \\Vu t \\ L 2 + l); 



(3. 
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/ 



p u u • Vu • u t dx 



p t u t ■ Vu • u t dx 



< \\p t t\\ L 2\\u ■ Vu|| L3 ||Ui|| L 6 
<C\\p tt \\ 2 L 2 + C\\Vu t \\ 2 L2 ; 
< l|pt||^K 2 |k3||Vu|| L6 < C||Vu t ||| 2 ; 



p t u ■ Vu t • u t dx 



< \\pt\\ L 3\\u\\ L oo\\Vu t \\ L 2\\u t \ 



IP 



< C\\Vu t \\v\\u t \\ L6 < C||Vu t ||2 2 . 

Substituting the above estimates back to (|3,32p and then integrating over [0, t] with 
< t < T, one has 



/ / puf t dxds + y sup / \Vu t \ 2 dx 
Jo Jr3 4 <t<T Jr3 

puf t dxds + C / / || Vu t \ 2 dxds + C. 
Jo Jm.3 



< Ce 



Here we used the following estimates: for all < t < T, 



(3.33) 



(p 7 )tdivu t dx 



< ll(^ 7 )t||L2||VUt|| L 2 



<C + ^ sup 

8 0<t<T 



|Vuj| dx; 



(pFF T ) t : Vu t dx 



< (\\Pt\\ L * + INb)||Vu t || i2 

V / lr -., ,2 



<C + ^ sup 

o 0<t<T 



|Vu(| dx; 



I 



p t \u t \ 2 dx 



div(pu)u^dx 



puVufdx 



and 



< \\^/pu t \\ L 2\\Vu t \\ L 2 < C + ^ sup / \Vu t \ 2 dx; 

o 0<t<TJR3 

/ /9 t (u • Vu)u((ix 

JR 3 



< ||/0t||x,s l|u|| X 6 ||Vll||x6 || Vll^ ||i2 



<C + ^ sup 

o 0<t<T 



\Vu t \ 2 dx. 



Choosing a sufficient small e, then (|3.23p . combining with GronwalFs inequality, gives 



I I pu 2 t dxdt + sup / \Vu t \ 2 dx < C. 

Jo JR 3 0<t<T Jr3 



The proof is complete. 

With aid of the above estimate, we can obtain the H 3 regularity as follows. 



□ 
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Lemma 3.8. 

IHUfffJ + ll F llL t °°//| + ll u llL t °°H3 < C. 
Proof. It follows from Lemma 13.61 and Lemma |3, 71 that 

GeLfHl Vp<eL?Hl div(pFF T ) G L t °°f£, 

which gives 

fiAu + (A + /x)Vdivu = G + AVp~f - div(pFF T ) G 
From the standard estimate for elliptic equations, we have 

||u|U°°H| < c. 

From the second equation in (jl.ip . Lemma 13, 61 and Lemma 13.71 we have 

pAu t + (A + p)Vdivu t = G t + AV(p^) t - div(pFF T ) t G L 2 ((0,T) x M 3 ). 

Therefore, again from the standard estimate for elliptic equations, we have 

u t G LfH%, G G £ 2 # 2 . 

With those bounds, we can apply a similar argument in Lemma 13.61 to deduce from the 
first and the third equation in (jl.ip that 

The proof is complete. □ 

Now we are able to state the proof of Theorem 13.11 

Proof of Theorem \3.1\ From the previous lemmas, the functions (p, F, u)\t=T* satisfy the 
same conditions as the initial data. Moreover 

G = pu t + pu • Vu G L^Hl 

and 

(-/iAu - (A + /x)Vdivu + AVp 7 )| i=T * = p(T*)g 

with 

^^(R 3 ), and y/p(T*)g G L 2 . 

Therefore, we can take (p, F, u)(T*) as the initial data and apply the local existence theo- 
rem to extend our local classical solution beyond T*, and this contradicts the assumption 
on T*. The proof is complete. □ 
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